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Abstract. The jets of microquasars with high-mass stellar companions are exposed to the dense matter field of the stellar 
wind. We present estimates of the gamma-ray emission expected from the jet-wind hadronic interaction and we discuss the 
detectability of the phenomenon at high energies. The proposed mechanism could explain some of the unidentified gamma-ray 
sources detected by EGRET instrument on the galactic plane. 
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1. Introduction 

Microquasars (MQs) are X-ray binary systems that present 
non-thermal jet-like features (Mirabel & Rodriguez 1999). The 
compact jets can be detected at radio wavelengths with flat 
spectra. The presence of apparent superluminal movements in 
some cases (e.g. Mirabel & Rodriguez 1994) indicates the exis- 
tence of relativistic bulk motions in the jet flow. Very recently, 
Chandra observations have revealed X-ray synchrotron emis- 
sion from the jet of XTE J 1550-564 (Corbel et al. 2002, Kaaret 
et al. 2003), a fact that indicates the presence of extremely rela- 
tivistic electrons with TeV energies and shock reacceleration in 
the jet. The interaction of such highly energetic particles with 
the ambient photon fields will produce high-energy inverse 
Compton (IC) radiation. In a series of recent papers, Aharonian 
& Atoyan (1998), Atoyan & Aharonian (1999), Markoff et 
al. (2001), Georganopoulos et al. (2002), and Romero et al. 
(2002) have studied different aspects of the synchrotron and 
the IC emission from MQs whose jets interact with different 
photon fields. Recently, Paredes et al. (2000) suggested that the 
MQ LS 5039 might be responsible for the gamma-ray source 
3EG J1824-1514. Kaufman Bernado et al. (2002) proposed 
that high-mass MQs with jets forming a small viewing angle 
(i.e. microblazars) might be the parent population of the set 
of bright and variable gamma-ray sources detected by EGRET 
along the galactic plane (e.g. Romero et al. 1999, Torres et al. 
2001). The main gamma-ray production mechanism for MQs 
advocated by Kaufman Bernado et al. (2002) was IC upscat- 
tering of UV photons from the stellar companion. However, 



pair creation absorption processes in the disk and coronal X- 
ray field might quench the source above a few MeV in many 
cases (Romero et al. 2002). 

In this Letter we present a new mechanism for the gen- 
eration of high-energy gamma-rays in MQs that is based on 
hadronic interactions occurring outside the coronal region. The 
gamma-ray emission arises from the decay of neutral pions 
created in the inelastic collisions between relativistic protons 
ejected by the compact object and the ions in the stellar wind. 
The only requisites for the model are a windy high-mass stellar 
companion and the presence of multi-TeV protons in the jet 1 . 
The presence of relativistic hadrons in MQ jets like those of SS 
433 has been inferred from iron X-ray line observations (e.g. 
Kotani et al. 1994, 1996; Migliari et al. 2002), although direct 
and clear evidence exists only for this source so far. In what 
follows we describe the model and present the results of our 
calculations. 



2. The jet 

The general situation discussed in this paper is shown in Fig.^ 
A binary system is formed by a black hole and a high-mass 
early-type star. A relativistic e - p jet is ejected by the compact 
object perpendicularly to the accretion disk plane. For simplic- 
ity, we shall assume that this is also the orbital plane, but this 
condition can be relaxed to allow, for instance, a precessional 
motion (Kaufman Bernado et al. 2002) or more general situa- 
tions as discussed by Maccarone (2002) and Butt et al. (2003). 
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1 Interactions of hadronic beams with moving clouds in the context 
of accreting pulsars have been previously discussed in the literature 
by Aharonian & Atoyan (1996). For an early discussion in a general 
context see Bednarek et al (1990). 
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Fig. 1. Sketch of the general situation discussed in the paper. A 
relativistic e - p jet is injected close to the black hole in a MQ 
with a high-mass stellar companion. The jet must traverse the 
matter field created by the stellar wind. The resulting interac- 
tion produces gamma-ray emission. In the figure perpendicu- 
larity is assumed between the jet and the orbital plane, but this 
particular assumption can be relaxed in a more general situa- 
tion. 

The jet axis, z, is assumed to be normal to the orbital radius 
a. We shall allow the jet to expand laterally, in such a way that 
its radius is given by R(z) = £z e , with e < 1 and Zo < Z < z m ax- 
For e = 1 we have a conical beam. The jet starts to expand at a 
height zo ~ a few hundred km above the black hole, outside the 
coronal region. The particle spectrum of the relativistic e — p 
flow is assumed to be a power law p (E' e ) = K e p E'~ a p , 
valid for E' e p min < E' e p < E' e p max , in the jet frame. The cor- 
responding particle flux will be J' e p (E' e p ) - (c/4-n)N' e p {E' e ^ p ). 
Since the jet expands, the proton flux can be written as: 



An K \ z ) 



F i-a 



(1) 



justifying the assumption that both directions are similar. Note 
that only photons emitted with angles similar to that of the in- 
clination angle of the jet will reach a distant observer, and thus 
8 can be approximated by the jet inclination angle. 

In order to determine the matter content of the jet we will 
adopt the jet-disk coupling hypothesis proposed by Falcke & 
Biermann (1995) and applied with success to AGNs (see also 
Falcke & Biermann 1996), i.e. the total jet power scales with 
the accreting rate as Qj = q^M&^c 1 , with qj = 10 -1 - 10 3 . 
The number density no' of particles flowing in the jet at Rq - 
R(zo) is then given by cnR^no' = Qj/m p c 2 , where m p is the 
proton rest mass. From here we can obtain no'. Additionally, 

n' = ClN'piE'p, zo) dE p '. Then, if £'™ x » E'f'\ which 

p 

is always the case, we have Ko = no {a - l)(E'™ m ) a ~ l , which 
gives the constant in the power-law spectrum at zo- 



where « > (a value n - 2 corresponds to the conservation 
of the number of particles, see Ghisellini et al. 1985), and a 
prime refers to the jet frame. Using relativistic invariants, it can 
be proven that the proton flux, in the observer (or lab) frame, 
becomes (e.g. Purmohammmad & Samimi 2001) 

T(F m cK 0(Zor ^(E p -p b jEl-rnyco S ey 
J p (E p ,0) = : 7777^> W 



sin 2 6» + r 2 cos6»- 



1/2 



where F is the jet Lorentz factor, 6 is the angle subtended by 
the emerging photon direction and the jet axis, and /3b is the 
corresponding velocity in units of c. The exponential depen- 
dence of the cross section on the transverse momentum (p t ) 
of the incident protons beams the gamma-ray emission into an 
angle <p < cp t /m p F ~ 0.17/F along the proton direction, hence 



3. The wind 

Early-type stars, like OB stars, lose a significant fraction of 
their masses through very strong supersonic winds. Typical 
mass loss rates and terminal wind velocities for O stars are 
of the order of 10~ 5 Mq yr~' and 2500 km s _1 , respectively 
(Lamers & Cassinelli 1999). At the base of the wind, the den- 
sity can easily reach 10 -12 g cm~ 3 . Such strong winds provide a 
field of matter dense enough as to produce significant hadronic 
gamma-rays when pervaded by a relativistic beam. 

The structure of the matter field will be determined essen- 
tially by the stellar mass loss rate and the continuity equation: 
M* = Anr 2 p{r)v(r), where p is the density of the wind and v 
is its velocity. The radial dependence of the wind velocity is 
given by (Lamers & Cassinelli 1999): 



v(r) = Voo ^1 - 



(3) 



where Vco is the terminal wind velocity, r» is the stellar radius, 
and the parameter /3 is ~ 1 for massive stars. Hence, using the 
fact that r 2 - z 2 + a 2 and assuming a gas dominated by protons, 
we get the particle density of the medium along the jet axis: 



n(z) 



4mn p Voa(z + a 2 ) 



+ a z 



(4) 



The stability of a relativistic jet under the effects of an ex- 
ternal wind has been recently investigated by Hardee & Hughes 
(2003) through both theoretical analysis and numerical sim- 
ulations. Their results indicate that jets surrounded by out- 
flowing winds are in general more dynamically stable than 
those surrounded by a stationary medium. Note that protons 
pertaining to the wind can diffuse into the jet medium. The 
wind penetration into the jet outflow depends on the parame- 
ter vj ~ vR{z)/D, where v is velocity of wind, R(z) is the radius 
of the jet at a height z above the compact object, and D is the 
diffusion coefficient, m measures the ratio between the diffu- 
sive and the convective timescale of the particles. In the Bohm 
limit, with typical magnetic fields Bo ~ 1-10 G, m < 1, and 
the wind matter penetrates the jet by diffusion. 
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4. Gamma-ray emission 



Table 1. Basic parameters of the model 



Pion decay chains leads to gamma-ray and neutrino production. 
The differential gamma-ray emissivity from 7r°-decays is: 



2Z ( 



(a) 



q y (E y ) = Ancr P p(E p )- 



a 



Jp(E 7 ,0)rj A . 



(5) 



Here, the parameter 77 a takes into account the contribution from 
different nuclei in the wind and in the jet (for standard compo- 
sition of cosmic rays and interstellar medium rj A = 1.4 - 1.5, 
Dermer 1986). J p {E y ) is the proton flux distribution evaluated 
at E — E y . The cross section o- pp (E p ) for inelastic p - p inter- 
actions at energy E p w 10E y can be represented above E p » 10 
GeVby cr pp (E p ) * 30x [0.95 +0.06 log(£ p /GeV)] mb. Finally, 
-*tP * s t ' le so " ca ^ ec ' spectrum-weighted moment of the inclu- 
sive cross-section. Its value for different spectral indices a is 
given, for instance, in Table Al of Drury et al. (1994). Notice 
that q y is expressed in ph s erg 1 when we adopt CGS units. 

The spectral gamma-ray intensity (photons per unit of time 
per unit of energy-band) is: 



Iy(Ey,9) 



f n(r' 
Jv 



)qy(r')d 3 r', 



(6) 



where V is the interaction volume. 

Since we are interested here in a general model and not in 
the study of a particular source, the spectral energy distribution 

-j 

Ly {E y , 9) = EyIy(E y , 6) is a more convenient quantity than the 
flux. Using eqs. Q, ©, (0 and l|6), we get: 

p^jfi a - 1 



Ly{E y ,9) 



2nm 1 p v tx 



(«7 



M,M disk o-pJIOEy) 



F- a+1 (E 7 -j3 h ^Ey-mjc 4 cos ej 



sin 2 6» + F 2 [cos 9- -M± 



1/2 



Vz 2 + a 2 



dz. 



(7) 



This expression gives approximately the 7r°-decay gamma-ray 
luminosity for a windy MQ at energies E y > 1 GeV, in a given 
direction 9 with respect to the jet axis. 

Depending on the characteristics of the primary star and the 
geometry of the system, high-energy gamma rays can be ab- 
sorbed in the anisotropic stellar photon field through pair pro- 
duction, and then inverse Compton emission from these pairs 
can initiate an e^pair cascade. This effect has been studied in 
detail by Bednarek (1997). The main effect of these cascades is 
a degradation of TeV gamma-rays into a form of softer MeV- 
GeV emission. Very close systems (a ~ 10 11 cm) with O stars 
and perpendicular jets can be optically thick for gamma-rays 
between ~ 0.1 - 100 TeV. If the jet is inclined towards the star, 
the effect can be stronger (Bednarek 1997). For systems with 
larger separations, the opacity rapidly falls below unity. 

In order to make some numerical estimates, we shall adopt 
the specific MQ model presented in Tabled The values chosen 
for the different parameters are typical for MQs with O stellar 
companions, like Cygnus X-l. We shall consider a conical jet 



Parameter 


Symbol 


Value 


Type of jet 


e 


1 


Black hole mass 




10 Mq 


Injection point 


Zl) 


50 Rl 

g 


Initial radius 


Ro 


5^, 


Radius of the companion star 


r* 


35 % 


Mass loss rate 


M 


10~ 5 Mq yr 1 


Terminal wind velocity 




2500 km s-' 


Black hole accretion rate 


M disk 


10~ 8 Mq yr 1 


Wind velocity index 


p 


1 


let's expansion index 


n 


2 


let's Lorentz factor 


r 


5 


Minimum proton energy 


rrf min 

E p 


10 GeV 


Maximum proton energy 


rrr max 

E p 


100 TeV 


Orbital radius 


a 




% = GM^/c 1 . 
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Fig. 2. Spectral high-energy distribution for windy MQs with 
proton index a = 2.2 and a = 2.8, for different jet/disk cou- 
pling constants (qj). The jet inclination with respect to the line 
of sight is assumed to be 10 degrees. An angle of 30 degrees 
reduces the luminosity in about two orders of magnitude. 



(e — 1) with conservation of the number of protons (n = 2) 
and a high-energy cutoff for the population of relativistic pro- 
tons of £ , '™ ax = 100 TeV. The minimum distance from the jet 
to the primary star is a — 70 ^0 ~ 5 x 10 12 cm. In Fig.|5]we 
show the spectral high-energy distribution for models with pro- 
ton index a = 2.2 and a = 2.8, and for different values of the 
jet/disk coupling parameter q ]7 with the results obtained using 
a numerical integration routine. We have added an exponential- 



like cutoff at E y 



0.1£7 ax . 



Since the forward momentum of 
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the protons in the jet is so great, the gamma-rays will be highly 
beamed, within an angle ~ arctg (R/z). Hence, only hadronic 
microblazars would be detected. 

5. Discussion 

Neutrinos are also generated by pion production chains. The 
signal-to-noise (S/N) ratio for the detection of such a v-signal 
can be obtained analyzing the event rate of atmospheric v- 
background and comparing it with the event rate from the 
source (see e.g. Anchordoqui et al. 2003 for details). The S/N 
ratio in a km-scale detector (like ICECUBE) in the 1-10 TeV 
band (including the effects of neutrino oscillations) is ~ 3 for 
one year of operation, assuming an inclination angle of 30 de- 
grees, a = 2.2 and that the neutrino spectrum roughly satisfies 
(e.g. Dar & Laor 1997): dF v /dE v 0.7dF y /dE y . This S/N is 
high enough as to justify speculations on the possibility of de- 
tecting a hadronic microblazar first from its neutrino signal (a 
serendipitious discovery in a detector like ICECUBE), and only 
later from its gamma-ray emission (through a pointed observa- 
tion). Note that this neutrino flux is different from, and in some 
cases can even be stronger than, the neutrino flux produced by 
photomeson processes within the coronal region (Levinson and 
Waxman 2001, Distefano et al. 2002), thus the detection of the 
neutrino emission from MQs seems promising. It is not possi- 
ble to separate in the neutrino signal the contributions form the 
photomeson and p — p channels, but simultaneous X-ray obser- 
vations can help to determine the characteristics of the relevant 
photon fields to which the inner jet is exposed, making then 
possible estimates of each contribution in particular cases. 

Another interesting aspect of the hadronic microblazar is 
that e ± are injected outside the coronal region through 7r ± de- 
cays. These leptons do not experience the severe IC losses 
that affect to primary electrons and pairs (Romero et al. 2002). 
These secondaries will mainly cool through synchrotron radi- 
ation (at X-rays, in the case of TeV particles) and IC interac- 
tions with the stellar seed photons (that would result into an 
additional source of MeV-GeV gamma-rays). The spectrum of 
secondary pairs roughly mimics the shape of the proton spec- 
trum. Hence, synchrotron emission from these particles will 
present indices a syn (a - l)/2, which for values of a ~ 2 
are similar to what is observed in MQs' jets at radio wave- 
lengths. The losses of the primaries in the inner source lead to 
a soft particle spectrum that is then injected in the region where 
the particles produce IC gamma-rays through interactions with 
stellar UV photons. Pure leptonic models for the gamma-ray 
flux (e.g. Kaufman-Bernado et al. 2002) then require particle 
re-acceleration in order to explain the flat radio-spectrum of 
sources like LS 5039 far from the core. In the hadronic model, 
the leptons are injected in situ with the right spectrum through 
hadronic decays, avoiding the problem. 

At TeV gamma-ray energies hadronic microblazars which 
are optically thin to pair production can be detected as uniden- 
tified, point-like sources with relatively hard spectra. This kind 
of sources could display variability. In the near future, new 
ground-based Cerenkov telescopes like HESS and MAGIC 
might detect the signatures of such sources on the galactic 
plane. Hadronic microblazars might be part of this population, 



as well as of the parent population of low latitude unidentified 
EGRET sources. 
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